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I. INTRODUCTION 



Among various theoretical attempts, the famous seesaw ideas provide us with a very 



natural way to understand why the masses of three 



mown neutrinos are so tiny compared to 



3y,Q,B,l6l,l7i. 



In the canonical type-I 



the masses of other Standard Model (SM) fermions 
seesaw mechanism, in which three right-handed (RH) neutrinos are introduced and lepton 
number violation is allowed, the effective mass matrix of three light Majorana neutrinos m^j 
is dramatically suppressed with respect to the electroweak scale if the RH neutrino mass 
matrix Mr is located not far away from the typical scale of grand unified theories. As 
a rough estimate, if light neutrino masses are stabilized around the sub-eV scale and the 
Dirac mass matrix between left- and right-handed neutrinos is comparable with the 
mass of the top quark, then Mr ~ 10^"^ GeV is naturally expected. The testability of such 
conventional seesaw models is therefore questionable. On the other hand, it is difficult to 
embed RH neutrinos into a theoretical framework at low-energy scales [e.g., the TeV scale to 
be explored by the Large Hadron Collider (LHC)] in a technically natural way while keeping 
the left-handed neutrinos to be light enough. For those realistically viable type-I or type-Ill 
seesaw models with the TeV-scale RH neutrinos, which could be experimentally accessible at 
the LHC, fine-tunings of cancellations among the contributions to rrii, from different heavy 
neutrinos have to be employed. These kinds of structural cancellations are usually attributed 
to some underlying flavor symmetries [s], o], Q, [ll|, Q. In the type-(I-l-II) seesaw model, 
one may also assume the mass term coming from RH neutrinos to be comparable with the 
one coming from the triplet Higgs, and thus, the light neutrino mass scale is brought down 



through a signiflcant cancellation between these two mass terms 



13|. However, to generate 



appreciable collider signatures of the heavy seesaw particles at the LHC, such a scheme 
potentially suffers from dangerous radiative corrections and requires unnatural flne-tuning 
even at loop level [l^ . 

In our recent work jl5], we have pointed out that the above-mentioned drawbacks of 
most TeV-scale type-I, type-Ill, and type-(I-l-II) seesaw models can be circumvented by 



considering the inverse seesaw model 



16|. In the latter framework, additional SM gauge 



singlets are adhibited together with a small Majorana mass insertion which explicitly breaks 
the lepton number. The phenomenology of this inverse seesaw mechanism is very rich: on 
the one hand, non-unitary neutrino mixing and CP violation can naturally show up and 
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are possible to be tested at the future long-baseline neutrino oscillation experiments; on the 
other hand, the heavy seesaw particles may result in very attractive signatures of lepton- 
flavor-violating (LFV) processes at the LHC Since the heavy singlets possess opposite CP 
signs and compose the pseudo-Dirac particles in pair, the lepton-number-violating (LNV) 
processes (such as the like-sign di-lepton events at the LHC) are significantly suppressed 
and practically invisible. 

However, the inverse seesaw scenario with three heavy singlet pairs contains too many 
free parameters and is not very predictive. To partly avoid this drawback, one may consider 
a simplified version of the generic inverse seesaw scenario by reducing its number of degrees 
of freedom. This sound motivation leads us to the minimal inverse seesaw scenario (MISS) 
at the TeV scale, in which only two pairs of the SM gauge-singlet neutrinos are introduced 
but the observed neutrino mass hierarchy and lepton flavor mixing can well be interpreted.^ 

The purpose of this work is to describe the MISS, which contains only two RH neutrinos 
(z/ri, Z/R2) and two SM gauge singlets {Si, S2), and to explore some of its low-energy con- 
sequences on neutrino mixing and CP violation. The remainder of our work is organized as 
follows. In Sec. HIl we will introduce the MISS and present some general formulas associated 
with the non-unitarity of the light neutrino mixing matrix. In Sec. Illlt we will constrain the 
parameter space of non-unitarity effects in the MISS by using current experimental data, 
calculate the Jarlskog invariants of leptonic CP violation, and discuss the discovery potential 
of such new CP-violating effects in the near detector of a neutrino factory. Finally, a brief 
summary will be given in Sec. llVi 



II. THE MINIMAL INVERSE SEESAW SCENARIO 

The MISS is constructed by extending the SM particle content with two RH neutrinos 
= i^Ri, ^R2) and two left-handed (LH) SM gauge singlets S = {Si, ^2). The mass part of 
the neutrino sector Lagrangian is then arranged so that it reads in the flavor basis 



-Cm= WMui^n + SM^ur + ^5^5^ + H.c. 



(1) 



where /i is a complex symmetric 2x2 matrix and Md and Mr are arbitrary 3x2 and 
2x2 matrices, respectively. Without loss of generality, one can always redefine the extra 



^ See also the minimal type-I seesaw ll|,ll8| and the minimal typc-(I+II) seesaw 19|. 
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singlet fields and work in a basis where fi is real and diagonal, namely, fi = diag(yU.i, yU.2) with 
fii < fi2- In addition to that, Mr can be made Hermitian by a further unitary transformation 
in the RH neutrino sector. The 7x7 neutrino mass matrix in the basis (i^l, i^^, S) is then 
rewritten as 



Md ^ 
ME Mi 
y Mr fi J 



(2) 



which is clearly a symmetric matrix with rank at most 6 20 1. 

Note that Mr is a SM singlet mass term, and hence, it is not governed by the scale of the 
SU{2)i^ symmetry breaking. In what follows, we will consider a particularly attractive case 
Mr > Md ^ yU, with Mr not far above the electroweak scale. In the limit fi ^ 0, the rank 
of Mj, reduces from 6 to 4, which leaves three light neutrinos massless. In reality, a tiny but 
non-vanishing /i can be viewed as a slight breaking of a global U{1) symmetry, and thus, it 
respects the naturalness criterion j21]. In contrast to the original inverse seesaw scenario, 
the MISS predicts one light neutrino to be exactly massless and brings in several interesting 
phenomena in neutrino oscillations and LFV processes. 

At the leading order in MdMj^\ the light neutrino mass matrix in the MISS is given by 



MbM^^h{MI)-^ME = FfiF^ 



(3) 



where F = MqMj^^ is a 3 x 2 matrix. For /i around the keV scale, MqMj^^ ~ 10~^ gives 
rise to the desired sub-eV light neutrino masses. In the limit /i — * 0, we have = 0, which 
corresponds to the lepton number symmetry restoration. The heavy sector consists of a pair 
of pseudo-Dirac neutrinos Pj = {Sj, i^rj) 2^ with a tiny mass splitting between the relevant 
CP-conjugated Majorana components of the order of fi. 

One can diagonalize mj, by means of a unitary transformation 



U'^rriyU* = rriy = diag(mi, 777,2, ms) 



(4) 



with rrii (for i = 1, 2, 3) denoting the light neutrino masses. We are left with either mi = 
(normal mass hierarchy Am^i > 0, NH) or 7773 = (inverted mass hierarchy Aml^ < 0, IH). 
Note that U itself is not the matrix that governs neutrino oscillations even if we choose a 
basis where the charged-lepton mass matrix is diagonal. 
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The LH neutrinos entering the charged- current interactions of the SM are superpositions 
of the seven mass eigenstates {umh, Pm) given at the leading order by 



where Ul^^M^^U^ = dia.g{m p^,mpj and iV ~ (l — |-^-^^) U 23|. Hence, one can write 



cc 
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-^W-i^r (iVz/„L + FUnPn^L) + H.c. 



(5) 



(6) 



The mixing between the doublet and singlet components in the charged currents results in 
several interesting phenomenological consequences: 

• The flavor and mass eigenstates of the left-handed neutrinos are connected by a non- 



unitary flavor mixing matrix N 



24( 1 ■ The magnitudes of non-unitarity effects in dif- 



ferent neutrino oscillation channels are predominated by the mass ratios between Md 
and Mr, and in principle, their underlying correlations have to be taken into account 
in analyses of the future experiments. 

The heavy singlets entering the charged currents due to the non-unitarity effects also 
enter the rare lepton decays, such as r yU7 and fi — > 67. Hence, unlike in the type-I 
seesaw model, their contributions to the LFV decays are not suppressed by the light 
neutrino masses 25|, but mainly constrained by the ratio M^M^^, which in principle 
admits observing these events in the upcoming LHC experiments. For example, the 
decays —>■ £^7 are mediated by P's and their branching ratios are given by 26| 



BR (4^^/37) 



2567r2M4.r„ 



i=l 



m 



Pi 



(7) 



where K = FUr, I (x) = -{2x^ + 5x^ - x)/[4(l - x)^] - 3xMnx/[2(l - x)^], and 
Ta is the total width of i^- In the conventional type-I seesaw model (i.e., without 
unnatural cancellations), one has approximately KK^ = 0{m^M^^), and therefore, 
BR {ia (^i3l) oc 0{m1) indicates a strong suppression of LFV decays. However, in 
the inverse seesaw model, one can have sizable K without any reference to the tiny 
neutrino masses, since the two issues are essentially decoupled. Thus, appreciable LFV 



rates could be obtained even for strictly massless hght neutrinos 



27| 



If the masses of the heavy singlets Pm do not fall far beyond the electroweak scale, in 
the MISS, as in the type-I seesaw model, one can expect an on-shell production of P 



at the LHC via the gauge boson exchange diagrams. The most distinctive signature 
would be the observation of LFV processes involving three charged-leptons in the final 
state 



In comparison to the generic inverse seesaw scenario, the MISS is certainly more restrictive 
and predictive because of its much fewer free parameters and underlying correlations among 
physical observables, which could be well tested at a future neutrino factory. In the remaining 
part of this work, we will discuss the phenomenological consequences mentioned above in 
more detail and, in particular, concentrate on the possible non-unitary CP-violating effects 
in neutrino oscillations. 



III. NON-UNITARITY EFFECTS 



A. Constraints on non-unitarity parameters 



As we have already mentioned, light neutrino mass eigenstates are connected to their 
flavor eigenstates by a non-unitary mixing matrix A^. In terms of the parametrization 



advocated in Ref. 29|], the non- unitary leptonic mixing is written as A = (1 — 77) f/ where 
the relevant Hermitian matrix 77 is given by r/ ^ ^^'^ the unitary matrix U can be 

parametrized in the standard form 



U 



C12C13 



■S12C13 



~Si2C23 — Ci2S23Si3e'^ C12C23 — Si2S23Si3e^'' 



Si3e 

■523^13 



■A 



■S12S23 ~ Ci2C23Si3e 



"C12S23 ~ ■Si2C23Si3e C23C13 



/1 



V 



V 



(8) 



provided Cij = cos% and Sij = sin^jj (for ij = 12, 13, 23). Since one light neutrino mass is 
vanishing in the MISS, we have only one Majorana phase. The current experimental bounds 
on 7] are rather stringent, implying that one can approximate ^jj's by the values of the 
mixing angles obtained from the neutrino oscillations. Present data on atmospheric, solar, 
and reactor neutrinos yield two neutrino mass-squared differences Amg^ — 7.65 x 10~^ eV and 



Amg^l ~ 2.40 X 10 eV, together with three neutrino mixings sin $12 ^ 0.304, sin^^: 



0.50, and sin^ ^ 



13 



0.01 



23 



30l |. There is no hint on the CP- violating phases, and thus. 



we take 6*13 < 10° and leave 5 as a free parameter in the following calculations. As for 
the non-unitarity parameters, |?7a^| are constrained mainly from universality tests of weak 
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interactions, rare leptonic decays, invisible width of the Z-boson and neutrino oscillation 
data. The present bounds on Irjapl (at the 90 % C.L.) are |31| : 



\V\ = ilVafsl) < 



2.0 X lO-^' 6.0 X 10-5 1.6 X 10"3 
8.0 X 10"^ 1.1 X 10-3 
~ ~ 2.7x10-3 



(9) 



In order to study the connections among physical parameters in the MISS, we adopt the 



parametrization of F from the work in Ref. 



321 ]. namely 



F = Uy/rfi^R^fi- 



(10) 



where 



R 



\ 



for the NH case, and 



R 





cos z — sin z 
sin z cos z 



^cos z — sin 



sm z cos z 




(11) 



(12) 



for the IH case. Here z = a + i/3 is an arbitrary complex number with both a and f3 being 
real. For sake of simplicity, in what follows we will use the parameters r = fii/fi2 and 
e = Am2i/|Am|]^| ~ 0.42. In the NH case, we obtain 



F=J'-^ 



m3 











\ 



-J^C23Sz C22.Cz 



+ Oie) 



(13) 



which yields 



1 ms rlc^P + 



2 /i2 



2r 





S23C23 
S23C23 C23 



(14) 
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Thus, the only non-neghgible off-diagonal entry is 77^,-, and we have the relation |?7^^| 



I Vf^T I ~ I VtT I 

obtain 



because, to a good approximation, 6*23 — 45°. Similarly, in the IH case, we 




nil 



B 



\ 



-j^C23X C23Y 



(15) 



which translates into 
/i 

1 mi 



V 



\ 



(16) 



-C23 {^AX* - BY*) S23 {}AX* - BY*) ^ 

clsCrW+m -S23C23(i|Xp+|Fp) , 

sls{}.\X\' + \Yn ) 

where A = ci2C^ + e'''si2S^, B = -ci2Sz + e^^si2Cz, X = S12C2 - e'''ci2Sz, and Y = S12S2 + 
e^^Ci2Cz- Hence, in both cases, the approximately maximal leptonic 23 mixing implies |?7^^| ~ 
\v^lT\ - \Vtt\ and |?7e^| ~ |?7er|. 

The allowed regions of non-unitarity parameters are illustrated in Figs. [1] and [2l In total, 
we have eleven parameters out of which seven (/ii, /i2, a, (3, p, 6, 9rs) are essentially free while 
the remaining four {612, 623, Amji, Amgj^) are quantities fixed by their best-fit values. In our 
numerical analysis, we randomly generate points in the seven-dimensional parameter space. 
The points that are within the 90 % CL. upper bounds on the non-unitarity parameters 
r]ai3 are plotted in the figures, where we do not resort to the approximate Eqs. f|T3|) - f|T6|) but 
rather utilize the precise Eq. (fTOj) . As a result, 10'' points build up each plot. Using Eq. Q, 
one can easily estimate that fi2 ^ should be ensured if there is no strong structural 
cancellations in F. Hence, in our numerical calculations, we set |Am3]^|/yU2 < 0.01 as a prior 
in order to avoid unnatural fine-tuning among model parameters. 

It is important to make clear that the decrease of the point-density towards smaller 
values of Irjapl is a mere numerical artifact, since these regions correspond to a physically 
not very interesting situation, and hence, it does not make sense to scan over such areas 
thoroughly. On the other hand, the density reduction observed in the opposite (i.e., growing 
|?7q,/}|) directions, c.f., the two upper plots in Fig. [H provides a true physical information and 
all the following statements are based on such kind of physically relevant features. 

In Fig. [H we present the allowed regions for the absolute values of the relevant non- 
unitarity parameters. The upper-left plot shows the correlation between |?7^^| and |?7/ir|j 
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10' 10' 10' 10"* 10' 10' 10' 10' 10' 10"* 10' 10' 



FIG. 1: Correlations among various parameters governing the non-unitarity effects in the MISS. 
Red points denote the normal neutrino mass hierarchy, while blue points correspond to the case of 
the inverted mass hierarchy. Generic experimental constraints are indicated by the green dashed 
lines (and, for simplicity, any would-be correlations in their determination have been neglected). 

which is stronger in the IH case than in the NH case. This is expected, since the analytical 
approximation in the IH case is better than that in the NH case. The upper-right plot shows 
the CP- violating phase of the parameter rj^r-, which is centered around in the NH case 
and around ±180° in the IH case, in agreement with Eqs. fll4l) and (1161) . Thus, in both 
cases, it is hard to achieve a sizable |?7^r| and maximal CP- violating effects simultaneously. 
The bounds on the other rj parameters are shown in the plots in the second row of Fig. [H 
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FIG. 2: Constraints on the underlying model parameters. As before, red and blue points denote 
the normal and inverted neutrino mass hierarchies, respectively. As described in the text, the 
\(3\ ^ 1 and /i2 — > regions correspond to a fine-tuning, which we wish to avoid, and thus, the 
scanning granularity is large. Moreover, ^ui ^ leads to a further reduction of the rank of rrii,, 
and hence, it is also disfavored. 

One can observe that \r]ee\ is mainly constrained in the NH case, while Irjrrl is restricted in 
the IH case. The allowed region for |?7eT-| is limited in both NH and IH cases, and there is no 
upgraded bound on |?7e^|. Note that, in the NH case, the leading order formulae f|T3|) and 
(IT^ do not provide as good approximation of the precise results as do the Eqs. (IT^ and 
fITB]) in the IH case, because the relevant effective expansion parameters (i.e., e"^ in the latter 
whilst only e in the former case) are not the same. This also imprints into the different 
widths of the red and blue bands in Fig. [H 

The allowed regions for /ij and z are shown in Fig. [2l There is no strong constraint on 
a. However, f] is bounded, since a larger (3 corresponds to a more sever fine-tuning of the 
model parameters. As for the yUj parameters, there are no generic upper bounds one should 
impose; however, smaller values of /ij correspond to a stronger fine-tuning entangled in F. 

Finally, let us summarize the upgraded bounds on the non-unitarity parameters we have 
obtained in the MISS under consideration: |?7ee| ^ 5.0 x 10~^, {rjerl ^ 8.9 x 10~^ in the NH 
case, and < 4.6 x 10"^ < 7.9 x 10-^ |r/^^| < 8.8 x 10"^ in the IH case. They 
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could be helpful when building a specific and realistic model based on the MISS. 



B. Jarlskog invariants 

In general, there are nine independent rephasing-invariant quantities that one can build 
at the quartic level out of the entries of a generic 3x3 lepton mixing matrix V {33|. 

4 = lm{Va^Vf,JV:^V^\) , (17) 

where the indices a ^ (3 run over e/i, /ir and re, while i ^ j can be 12, 23 and 31. 
Note that all nine J*{'s coincide if is a unitary matrix, since all six unitarity triangles, 



despite their different shapes, span the same area [3J, |35|]. However, if non- unitarity effects 
(like those studied in this work, i.e., V = N) are present, this is no longer the case and 
one can expect deviations from such a simple picture driven by the relevant non-unitarity 
parameters (denoted by 77^^ in the current study). In such a case, it is instructive to know 
which configuration oi a ^ /3 and i ^ j is most affected for a specific non-unitarity pattern. 
In particular, to leading order in rjajj, one can write 

4^J + A4, (18) 

where J = ci2C^3C23Si2Si3S23 sin 5 governs the CP- violating effects in the unitary limit. The 
second term in Eq. (fT8|) depends on the off-diagonal (generally complex) r/'s, and hence, it 
does not necessarily vanish in the limit Oi^ and might even dominate the CP-violating 
effects. The complete expressions for AJ^^^^ are listed in Appendix \M Focusing on the 
dominant off-diagonal entry |?7^r|; c.f. Fig. [H the following two contributions survive for 
?7e^ and 6*13 0: 

AJ^l = -|?7^^| sin5^^sin26'23COS^6'l2 - |'7er|Sl2Cl2S23C23Sin(5er , (19) 

AJ^l = |?7^^| sin5^^sin26'23sin^6'i2 - |v7er|si2Ci2S23C23sin5e^ , (20) 

provided 77^^ = Irja/sle^^"'^ ■ Even beyond the simple limit above, one can observe another 
interesting feature: 

t23 _ t31 _ t23 _ t31 /9-1 \ 

where all the small parameters {rjaiB and S13) have been kept at linear order. However, these 
relations are mere reflections of the smallness of the mixing angle 6*13. 
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FIG. 3: Constraints on the Jarlskog invariants. Here, only the NH case is considered, since the 
CP-violating effects due to the non-unitarity parameters are small in the IH case. It is so partly 
due to the smallness of rjer and also because of the preferred value of the 5^r phase, c.f. Fig. [TJ 



In Fig. [3], we illustrate the correlations between Jarlskog invariants in the NH case. As 
we expect, they are linearly dependent, and the spread (i.e., the finite width of the allowed 
strips) is due to higher-order corrections. Although the magnitude of AJ^l does not seem to 
be comparable to that of J, we will show later that the CP-violating effects induced by the 
phases of non-unitarity parameters can be quite significant, since they are not suppressed 
by the small neutrino mass-squared difference Am'^i- For the IH case, according to Eq. ( |T6l) . 
Vefj. ~ Ver are small quantities, and the phase of r/^r is close to vr, c.f. Fig. [H Hence, there is 
no observable CP-violating effect coming from the non-unitarity parameters in the IH case. 



C. Sensitivity search at a neutrino factory 



For a non-unitary lepton flavor mixing inatrix N, the vacuum neutrino oscillation tran- 
sition probability Pq,^ can be written as 36| 



i,3 



i>j 



t>j 



2E 



12 



where Am% 



are the neutrino mass-squared differences and JF* are defined by 



7.P 



(23) 



with the normahzed non-unitary factor 

Rap ~ 



v) 



al3 



(24) 



If Earth matter effects are taken into account, then one can replace the vacuum quantities 
U and rrii by their effective matter counterparts, see e.g . Ref. 37|. 



As mentioned in the hterature 



38 



39 



40 



41 



42 



431], the i^fj, —>■ Vr channel together with 



a near detector located at a short distance provides the most favorable setup to constrain 
the non-unitarity effects.^ In this respect, we consider the transition probability P^r for 
a neutrino factory with a sufficiently short baseline length L. We neglect the tiny matter 
effects and small contributions of 6*13 and Am|^. Then, P^^ reads 



^ p,T 



4s^3C^3Sin2 ' ^"^31^ 



A\ri^^\ sm(3^^S23C23 sm ' 



(25) 



AE J ^' \ 2E 

where E is the neutrino beam energy and the second term is CP-odd due to t he phas e (5^,-, 



and hence, distinctive CP-violating effects can appear in neutrino oscillations [29|, |46|] . The 
last term in Eq. fl25l) plays the dominant role at 'zero' distance, as it does not depend on L. 

Note that the shape of the CP-odd term is justified by the structure of the relevant 
Jarlskog invariants derived in Sec. IIII B[ Indeed, the Im(jF^gjF^o) factors in Eq 



cor- 



respond to — [apart from the irrelevant real rescaling of J due to the denominator in 
Eq. (1241) ]. This means that the sum over i > j in the last term of Eq. (!22|) is proportional 
to J^l — J"^ = \r]fj_r \ sin sin 26*23 provided Amg^ ^ —Aml^. 

In order to show the feasibility of observing such a signal in the future long-baseline 
neutrino oscillation experiments, we consider a typical neutrino factory setting with an 
OPERA- like near detector with fiducial mass of 5 kt. We assume a setup with approximately 
10^^ useful muon decays and five years of neutrino running and another five years of anti- 



neutrino running. We make use of the GLoBES package 



47 



48l ] with a slight modification 



of the template Abstract Experiment Definition Language (AEDL) file for the neutrino 



An alternative study for the disappearance channel Vfj, —f together with a far detector located at 



7500 km has been performed in Refs. 4Jjl45|- In particular, a sensitivity of 0(10 ^) could be achieved 



due to matter effects, for which the neutrino oscillation probability is only linearly suppressed in r/^^. 
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FIG. 4: Sensitivity limits at 90 % C.L. on the non-unitarity parameter r/^,- as a function of the 
basehne length L. Solid curves denote the parent muon energy = 50 GeV with CP phases 
being labeled in the figure, while the dashed curve corresponds to i?^ = 25 GeV and 5^r = 0. The 
bounds on |?7^t| in the MISS for the NH and IH cases are also shown by red lines. 



factory experiments [49|, |50|. In Fig. HI we display the sensitivity to function of 

the baseline length L for the near detector. One can observe that such a setup provides 
indeed an excellent probe for this type of non-unitarity effects. An interesting feature, 
which appears if 5^^ is sizeable, is that the sensitivity of the near detector will be improved 
around the baseline length L ~ 300 km. This sensitivity enhancement could mainly be 
regarded as a compromise between new physics effects and the standard neutrino oscillation 
behavior. At a very short distance, the transition probability P^r is determined by the 
last term of Eq. (l25l) . whereas with increasing L, the second term gradually dominates the 
flavor transitions. Thus, a distance L < 500 km (i.e., the CERN-Frejus distance) would be 
favorable for the near detector. 



IV. SUMMARY 

We have proposed the MISS — an economical low-scale seesaw scenario with minimal 
particle content in the framework of the inverse seesaw model. Compared to the generic 
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inverse seesaw mechanism, only two pairs of SM gauge singlets are introduced into the MISS, 
which gives rise to strong correlations among the non-unitarity parameters. Since one light 
neutrino has to be massless in this scenario, we have discussed the experimental constraints 
on these non-unitarity parameters in both NH and IH cases. In view of our numerical and 
analytical results, the only possibly sizable and phenomenologically interesting non-unitarity 
parameter is r/^r, and the current upper bound on |?7^t-| is improved from 1.1 x 10"'^ to 
7.9 X 10~^ in the IH case. The Jarlskog invariants in the presence of non-unitary neutrino 
mixing have been calculated. The relative CP-violating phase of 7]^^ is well constrained 
by the structure of the MISS, and there are essentially no observable CP-violating effects 
induced by S^r in the IH case. We have also shown that the CP-violating effects emerging in 
the MISS can be well tested at a future neutrino factory with an OPERA-like near detector 
at a distance less than a few hundred kilometers. The possible collider signatures at the 
LHC and the LFV processes are also promising. However, a detailed analysis exceeds the 
scope of the current work, and hence, it will be elaborated on elsewhere. 

Finally, we would like to stress that the MISS is motivated not only by its simplicity and 
predictivity, but it can in particular be viewed as a limit of the "standard" ISS setting when 
the heaviest pseudo-Dirac neutrino essentially decouples. This is the kind of behavior one 
would expect in grand unified models where a hierarchy in the RH sector is often natural 
because of the link of the RH-triplet Yukawa couplings to the other parts of the Yukawa 
sector.'^ Alternatively, it is often realized in certain classes of flavor symmetric models in 
which the RH neutrinos are assigned to some two-dimensional representations of the flavor 
group, i.e., the smallest group containing one-, two-, and three-dimensional representations 



of the symmetric permutation group 5*4 53|]. One can also adopt a variant of the several 



strategies proposed to govern the flavor structure of the original inverse seesaw framework. 



in particular, to accommodate the tri-bimaximal mixing pattern, see e.g. Ref. [5J] and 
references therein. However, a thorough implementation of a favor symmetry in the given 
context is beyond the scope of the current work, and will not be further discussed here. 



3 



In this respect, it is worth noting that the smaUness of the RH neutrino mass scale does not in general 



obstacle the grand-unified constructions, see for instance Refs. 



51 



and references therein. 
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APPENDIX A: CALCULATION OF JARLSKOG INVARIANTS 

The Jarlskog invariants for a non-unitary lepton flavor mixing matrix are defined by 

J% = lm{N^,Np,K^N;,) , (Al) 

provided = {l—rfjU, where U and rj are 3x3 unitary and Hermitian matrices, respectively. 
The stringent experimental constraints on the deviation of A^ from U allow one to perform 
an expansion of Eq. (]A1|) in powers of the small parameters rjap and ^13. Up to the second 
order in 77^/3 and S13, one has J^-'j ^ J + AJ^-^^ , where J = ci2C^3C23Si2Si3S23 sin 5 and 



7 

+ v:,u^.Up,u;^u;, + v;,u^,Up,u:^u;,) . (A2) 

In the parametrization ([8]), the nine relevant Jarlskog invariants read: 

AJg^ = -|?7eM|Sl2Cl2C23(l + C23)sin(5e^-h |?7er|Sl2Cl2S23C23Sin5er , (A3) 

A^e^ = |?7e^|Sl2Ci2S23C23 sin 5e/x + |?7er|Sl2Cl2S23C23 singer , (A4) 

AJ^l = -AJ^f^-2\r]^r\cl2S23C23Sm6f,r , (A5) 

AJ31 = -Aj23 + 2|r/^,|4s23C23sin5^, , (A6) 

AJ^g = |77eMkl2Cl2S23C23 sin 5eAt - |?7er I Sl2Cl2S23(l + S23) singer , (AT) 

AJ2 = AJJ = -AJ- = AJ^^l = A41 . (A8) 
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